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Purpose: We examined expression profiles of 16 micro RNAs 
(miRNAs) in triple negative breast cancers to identify their poten- 
tial as biomarkers for lymph node metastasis. Methods: The ex- 
pression profiles of miR-9, miR-21, miR-30a, miR-30d, miR-31, 
miR-34a, miR-34c, miR-100, miR-122, miR-125b, miR-146a, 
miR-146b, miR-155, miR-181a, miR-200c, and miR-205 were 
examined by using real-time quantitative reverse transcription 
polymerase chain reaction in tumor samples and corresponding 
benign breast tissues. Their associations with histopathological 
features and prognostic parameters were assessed. Results: 
When compared with the expression in benign breast tissues, 
seven of the miRNAs (miR-31, miR-205, miR-34a, miR-146a, 



miR-125b, miR-34c, and miR-181a) were downregulated more 
than 1 .5-fold in tumor tissues, whereas, only miR-21 was found 
to be upregulated more than 1 .5-fold in tumor tissues. Although 
miR-200c levels were decreased only 1 .12-fold in tumor tissues, 
the reduced expressions of miR-200c and miR-205 were signifi- 
cantly associated with lymph node metastasis (p= 0.021 and p= 
0.016, respectively). Conclusion: Our results demonstrate that 
miR-205 and miR-200c expression levels may be useful in pre- 
dicting lymph node metastasis in triple negative breast cancer 
patients. 
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INTRODUCTION 

MicroRNAs (miRNAs) are small noncoding RNAs that 
regulate gene expression through base-pairing with their tar- 
get messenger RNAs. Despite the fact that their functions are 
largely unknown, most miRNAs have the potential to regulate 
more than one target gene in the human genome. Studies 
have revealed that miRNAs may become deregulated in hu- 
man cancers, acting as oncogenes or tumor suppressor genes 
[1,2]. miRNAs also regulate genes associated with prognostic 
parameters, which play a role in tumor progression, invasion, 
or metastasis. Oncogenic miRNAs (oncomirs) are usually 
overexpressed in cancer tissues, whereas tumor suppressor 
miRNAs usually demonstrate reduced expression [3]. 
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Breast cancer may be divided into four therapeutic and 
prognostic subgroups according to expression of the estrogen 
receptor (ER) and/or progesterone receptor (PR) and over- 
expression of human epidermal growth factor receptor 2 
(HER2) [4]. The ER and PR are expressed in the vast majority 
of tumors, and they are usually associated with a favorable 
prognosis. Overexpression of HER2 is associated with more 
aggressive behavior, but patients in this group do have the op- 
portunity to benefit from treatment with trastuzumab. In con- 
trast, patients in the triple negative breast cancer (TNBC) 
group have the poorest outcome, and these patients do not 
benefit from either endocrine treatment or trastuzumab. In 
our patient series, 50% of the TNBC patients were found to 
have lymph node metastasis at the time of surgery, a figure far 
higher than that of the other types, and one that reflects the 
aggressive behavior of this type of breast cancer. Although it 
represents a small fraction of breast cancer, TNBC is com- 
posed of a heterogeneous disease group with different histo- 
logic, biologic, and behavioral features [4-6]. 

In the present study, we aimed to investigate the expression 
profiles of 16 miRNAs for their utility as predictive markers of 
lymph node metastasis in TNBC. We also assessed their asso- 
ciation with clinical and histopathological parameters. 
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METHODS 

Patients 

Formalin-fixed, paraffin-embedded pretreatment tissue 
samples from 32 TNBC patients who underwent radical mas- 
tectomy and axillary dissection between the 2004 and 2011 
were included in the current study The study was conducted 
in accordance with the Declaration of Helsinki and approved 
by our Institutional Ethics Committee (IRB approval number, 
2013-107). None of the patients had a known history of famil- 
ial breast cancer. Baseline characteristics of the patients are 
presented in Table 1. The ER, PR, and HER2 status of tumor 
samples that were determined by immunohistochemical 
methods in routine pathological examination and the histo- 
pathological features of tumor tissues, were re-evaluated by 1 
pathologist in each case. In 12 of the cases, HER2 negativity 
was also confirmed by the fluorescence in situ hybridization 
(FISH) method. 

Selection of miRNAs 

The miRNA selections were made based on a literature 
search for previous reports. To identify miRNAs with a predic- 
tive potential for lymph node metastasis, the search was fo- 
cused on miRNAs which were found to 1) be related to the tri- 
ple negative phenotype; 2) play critical roles in the cell cycle; 3) 
have prognostic or predictive values; 4) have a role in the inva- 
sion-metastasis cascade; or 5) be associated with aggressive 
behavior in several malignancies including breast cancer. Out 
of several candidates, further selection was made according to 
the target genes of the miRNAs for their involvement in the 
tumorigenic pathways of breast cancer. In an attempt to focus 
on relatively less-studied miRNAs, we excluded miRNAs such 
as miR-lOb, miR-373, and miR-520c, the roles of which have 



Table 1 . Baseline characteristics of triple negative breast cancer (TNBC) 
patients 



Characteristic 


No. (%) 


No. of TNBC patients 


32 


Age (yr)* 


59.3±15.6 


Histology 




Ductal 


23 (71.9) 


Lobular 


4(12.5) 


Other 


5(15.6) 


Tumor size (cm) 




<2 


8 (25.0) 


>2 


24 (75.0) 


Lymph node metastasis 




Absent 


16(50.0) 


Present 


16(50.0) 


*Mean±SD. 
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been confirmed in the metastatic pathways of breast cancer [3]. 
As a result, out of hundreds of known human miRNAs, 16 
candidates were considered to be suitable for the study. Select- 
ed miRNAs and their primers are listed in Table 2. 

Isolation and relative quantification of miRNAs 

For miRNA extraction from cancer tissues, previously pre- 
pared tissue macroarray blocks containing only tumor areas 
were used. In all cases, corresponding normal breast tissues 
were also selected for miRNA extraction. 

Total RNA was isolated using a miRNeasy FFPE Kit (QIA- 
GEN, Hilden, Germany). The quantity of the RNA was as- 
sessed by using a NanoDrop 1000 spectrophotometer (Ther- 
mo Scientific, Wilmington, USA). miRNAs were amplified 
using the SYBR Green method after reverse transcription and 
cDNA synthesis from isolated total RNA. In these steps, miS- 
cript technology kits were used according to the manufactur- 
ers instructions (QIAGEN). RNU6B small nuclear RNA was 
used as a reference with serial dilutions of template cDNA. All 
reactions were performed in a final volume of 20 [iL contain- 
ing 2 |iL of synthesized cDNA and carried out on RotorGene 
Q-5 Plex (QIAGEN) in duplicate by including no template 
and reverse transcriptase minus controls. After the amplifica- 
tion steps, melting curves were generated from melting data 
acquired from 50°C to 85°C, at a ramping rate of FC/sec, to 
evaluate the polymerase chain reaction products. miRNA ex- 
pression levels were calculated using the AACt method [7]. 

Statistical analysis 

Results are presented as mean ± standard deviation for con- 
tinuous variables, and defined as percentages for categorical 



Table 2. The list of micro RNA (miRNA) and mature miRNA primers 
used in the study 



miRNA 




Mature miRNA primer 


Hsa-miR- 


■9-3p 


ATAAAGCTAGATAACCGAAAGT 


Hsa-miR- 


■21 -5p 


TAGCTTATCAGACTGATGTTGA 


Hsa-miR- 


■30a-5p 


TGTAAACATCCTCGACTGGAAG 


Hsa-miR- 


■30d-5p 


TGTAAACATCCCCGACTGGAAG 


Hsa-miR- 


■31 -5p 


AGGCAAGATGCTGGCATAGCT 


Hsa-miR- 


■34a-5p 


TGGCAGTGTCTTAGCTGGTTGTT 


Hsa-miR- 


■34c-5p 


AGGCAGTGTAGTTAGCTGATTGC 


Hsa-miR- 


■100-5p 


AACCCGTAGATCCGAACTTGTG 


Hsa-miR- 


■122-5p 


TGGAGTGTGACAATGGTGTTTG 


Hsa-miR- 


■125b-5p 


TCCCTGAGACCCTAACTTGTGA 


Hsa-miR- 


■146a-5p 


TGAGAACTGAATTCCATGGGTT 


Hsa-miR- 


■146b-5p 


TGAGAACTGAATTCCATAGGCT 


Hsa-miR- 


■155-5p 


TTAATGCTAATCGTGATAGGGT 


Hsa-miR- 


■181a-5p 


AACATTCAACGCTGTCGGTGAGT 


Hsa-miR- 


■200c-3p 


TAATACTGCCGGGTAATGATGGA 


Hsa-miR- 


■205-5p 


TCCTTCATTCCACCGGAGTCTG 
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variables. We used the Kolmogorov-Smirnov test to assess the 
distributions of continuous variables. Mann-Whitney non- 
parametric tests were used to compare the differences in ex- 
pression levels of miRNAs between different clinicopathologic 
disease groups such as metastatic and localized tumor groups. 
Unsupervised hierarchical cluster analysis was also performed 
to evaluate miRNA expression signatures in different tumor 
subtypes. All analyses were performed using SPSS for Win- 
dows software package, version 11.5 (SPSS Inc., Chicago, 
USA). All p-values were two-tailed, and p < 0.05 was set as the 
level of statistical significance. 

RESULTS 

To identify their predictive potential for lymph node metas- 
tasis, we examined the expression profiles of 16 miRNAs in 32 
TNBC samples and corresponding normal breast tissues. Sev- 
en of the miRNAs were downregulated in tumor tissues by 1.5- 
fold or more, and only miR-21 was found to be upregulated by 
more than 1.5-fold. The mean miRNA expression levels of 
cancer tissues normalized to normal breast tissues and the ra- 
tios in patients with and without metastasis are presented in 
Figure 1. miR-31 was the most strongly downregulated miR- 
NA (mean, 0.247 ± 0.299), followed by miR-205 (mean, 0.367 ± 




0.1 L 

i- LO 

± 1 
E 1 

Figure 1. Mean micro RNA (miRNA) expression levels of overall, metastatic and nonmetastatic triple negative breast cancer samples. Data were nor- 
malized to benign breast tissue levels and standardized by log2 transformation. Vertical lines indicate standard error. 
LN+ = lymph node metastasis; LN- = lymph node nonmetastasis. 
*p<0.05. 



0.353) and miR-34a (mean, 0.396 ±0.303), respectively (Table 
3). When the miRNA levels of the metastatic and nonmeta- 
static groups were compared, miR-205 and miR-200c levels 
were significantly decreased in the metastatic group (p- 0.016 
and p- 0.021, respectively). In the metastatic group, miR-205 

Table 3. Comparison of micro RNA (miRNA) expression levels in meta- 
static and nonmetastatic triple negative breast cancer 



Mean miRNA expression change (fold) 

miRNA p-value 

Metastatic tumors Nonmetastatic tumors 



Hsa-miR-9 


0.638 


0.821 


0.309 


Hsa-miR-21 


1.873 


1.691 


1.000 


Hsa-miR-30a 


0.753 


0.890 


0.274 


Hsa-miR-30d 


0.706 


0.869 


0.291 


Hsa-miR-31 


0.213 


0.280 


0.509 


Hsa-miR-34a 


0.385 


0.407 


0.546 


Hsa-miR-34c 


0.600 


0.682 


0.474 


Hsa-miR-100 


0.640 


0.751 


0.440 


Hsa-miR-122 


0.648 


1.260 


0.051 


Hsa-miR-125b 


0.547 


0.677 


0.546 


Hsa-miR-146a 


0.492 


0.523 


0.940 


Hsa-miR-146b 


0.728 


1.223 


0.228 


Hsa-miR-155 


1.319 


1.589 


0.910 


Hsa-miR-181a 


0.649 


0.738 


0.734 


Hsa-miR-200c 


0.591 


1.202 


0.021 


Hsa-miR-205 


0.265 


0.470 


0.016 
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Table 4. Predictive values of miR-205 and miR-200c expression chang- 
es using 5-fold decrease as treshold 





Predicted 


Predicted 


Sensitivity 


Specificity 




metastatic 


nonmetastatic 


No. (%) 


No. (%) 


miR-205 levels 










Observed 


11 


5 


11/16(68.8) 


13/16(81.3) 


metastatic 










Observed 


3 


13 






nonmetastatic 










miR-200c levels 










Observed 


3 


13 


3/16(18.8) 


16/16(100) 


metastatic 










Observed 


0 


16 






nonmetastatic 










miR-205+miR-200c 










levels 










Observed 


14 


2 


14/16(87.5) 


13/16(81.3) 


metastatic 










Observed 


3 


13 






nonmetastatic 











expression levels showed a more than 5 -fold decrease in 1 1 
(68.8%) cases, whereas this was observed in only three (18.8%) 
cases in the nonmetastatic group. miR-200c levels also exhibit- 
ed a more than 5-fold decrease in four (25.0%) of the metastat- 
ic cases, whereas none of the cases in the nonmetastatic group 
demonstrated such a decrease. Collectively, a 5 -fold decrease 
in either miR-205 or miR-200c expression levels was observed 
in 14 of 16 metastatic cases (87.5%). Using miR-205 and miR- 
200c expression levels, 22 of 32 patients (71.9%) could be cor- 
rectly classified according to lymph node metastatic status. Us- 
ing a 5 -fold decrease as the threshold for either of these miR- 
NAs, lymph node metastasis could be predicted with 87.5% 
sensitivity and 81.3% specificity (Table 4). Notably, miR-200c 
levels were decreased in the metastatic group (mean, 0.591 ± 
0.531), whereas these levels were unchanged in the nonmeta- 
static group (mean, 1.202 ± 1.009). Similar results were also 
observed for miR-122 and miR-146b, which exhibited more 
than 1.5 -fold higher levels in nonmetastatic tumors than in 
metastatic tumors; however, the differences between both 
groups did not reach the level of significance {p- 0.051 and 
p- 0.228, respectively). Except for miR-200c, miR-122, and 
miR-146b, all other miRNAs were either decreased or in- 
creased simultaneously in both the metastatic and nonmeta- 
static groups (Figure 1). miR-21 and miR-1 55 were the most 
highly upregulated miRNAs in both groups. 

The normalized ratios of tumor/benign breast tissue miRNA 
expression levels in each case are shown by heatmap colors in 
Figure 2. We did not observe a separation of metastatic and 
nonmetastatic tumors in hierarchical clustering of tumor sam- 
ples. However, miRNA expression levels revealed the presence 
of different tumor clusters, and most miRNA levels were high- 
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Figure 2. Unsupervised hierarchical clustering analysis result according 
to all micro RNA (miRNA) levels. Heatmap colors represent normalized 
miRNA levels in each case, and each line represents a case. Red boxes 
in the first column indicates lymph node metastasis (LNM). Data are 
standardized by log 10 transformation. 

ly correlated within these clusters. 

We did not observe any significant association between 
miRNA expression levels and histological type, histological 
grade, nuclear grade, or tumor diameter (data not shown). 
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DISCUSSION 

The accumulation of data from previous reports revealed 
that deregulation of miRNA expression may play significant 
roles in the progression and dissemination of breast cancers. 
In order to identify their predictive potential for lymph node 
metastasis, we selected candidate miRNAs that play critical 
roles in the cell cycle or in tumor initiation or progression. 

We found miR-31 to be the most strongly downregulated 
miRNA in TNBC tissues. miR-31 was previously shown to in- 
hibit multiple steps of the invasion-metastasis cascade in breast 
cancer [8]. Its expression is reported to be highly decreased in 
metastatic breast cancer cell lines, and decreased levels have 
also been observed in nonmetastatic breast cancer cell lines [3]. 
The promoter-associated CpG island of the LOC554202 non- 
coding gene, the host gene of miR-31, was found to be heavily 
methylated in TNBC cell lines, which may lead to the loss of 
miR-31 expression and the activation of its target prometasta- 
tic genes [9]. Although the difference was not significant, met- 
astatic tumors exhibited lower expression levels of miR-31 
than did nonmetastatic tumors in our study. These findings 
were in accordance with previous reports and support the tu- 
mor suppressor activity and inhibitory effects of miR-31 on 
the invasion-metastasis cascade in TNBC. 

miR-205 regulates the epithelial-mesenchymal transition 
(EMT), which is related to aggressive behavior, and it has been 
found to suppress the cell growth and invasive capacity of 
breast cancer cells [10,11]. miR-205 has been reported to be 
underexpressed in TNBC tissues in previous studies, a finding 
which supports its tumor suppressor function in breast cancer 
[12,13]. In our study, when compared with benign breast tis- 
sues, miR-205 expression levels were found to be decreased 
more than 2-fold in tumor tissues. Additionally, the decrease 
in the metastatic group was significantly higher than that in 
the nonmetastatic group in TNBCs. 

In the present study, the second significantly decreased miR- 
NA in the metastatic group was miR-200c. In contrast to miR- 
205, miR-200c levels were downregulated in metastatic tumors 
but the levels were not changed in nonmetastatic TNBCs. The 
MiR-200 family, including miR-200c, has been shown to in- 
duce the epithelial phenotype and to suppress EMT, as does 
miR-205. MiR-200c has been reported to suppress the tumori- 
genicity of breast cancer stem cells in vivo [14], and overexpres- 
sion of miR-200s reduced lung metastases in an experimental 
model of pulmonary metastasis, findings that suggest the cru- 
cial role of these miRNAs in the metastasis of breast cancer [3]. 
Downregulation of miR-200c was also found to be associated 
with the triple negative phenotype in breast cancers [15]. 

Because miR-205 and miR-200c levels were significantly de- 



creased in metastatic TNBCs, and their expression levels ex- 
hibited a trend with aggressive behavior, we stratified patients 
according to miR-205 and miR-200c expression levels using a 
5-fold decrease as a threshold. In this manner, our data re- 
vealed that lymph node metastasis could be predicted with 
87.5% sensitivity and 81.3% specificity in TNBCs. Although 
miR-205 was found to be more sensitive than miR-200c in 
predicting lymph node metastasis, it is worth noting that more 
decreased levels of miR-200c were only observed in metastatic 
tumors, and its levels were not changed in nonmetastatic tu- 
mors. Therefore, the results obtained with miR-200c in our 
study are promising for other members of the miR-200 family 
to be predictive markers of lymph node metastasis in TNBCs. 

In our study, miR-34a, miR-34c, miR-146a, and miR-181a 
levels were found to be decreased more than 1.5-fold in TNBCs. 
However, we did not find any significant differences in their 
expression levels between metastatic and nonmetastatic tu- 
mors. On the other hand, miR-146b and miR-122 levels were 
decreased in the metastatic group, whereas they were un- 
changed in the nonmetastatic group. Although the differences 
in their expression levels between the metastatic and nonmet- 
astatic groups did not reach the level of significance, the small 
number of patients limits the statistical power of our study, 
and thus the predictive values of miR-146b and miR-122 may 
be underestimated for this reason. 

miR-21 and miR-1 55 were the only upregulated miRNAs in 
the present study. Circulating levels of miR-1 55 have been re- 
ported to be associated with overt metastases in breast cancer 
patients [16]. However, in another study, miR-1 55 levels did 
not correlate with metastatic relapse in stage II and III patients 
[15]. We did not observe a difference between the miR-1 55 
levels of tumor tissues in metastatic and nonmetastatic pa- 
tients. Moreover, the increase in levels in tumor tissues was 
only 1.45 -fold when compared to that in benign breast tissues. 
MiR-21 has been found to be overexpressed in several human 
tumors, and its levels were also reported to be increased in 
TNBCs [12,17]. Our data revealed 1.8-fold overexpression in 
miR-21 levels in tumor tissues versus benign breast tissues. 
However, its levels were increased in both metastatic and non- 
metastatic tumors, and we did not observe a difference in ex- 
pression levels between the two groups. 

We did not observe any change in expression levels of miR- 
9, miR-100, miR-30a, and miR-30d in tumor tissues. miR-30a 
and miR-30d belong to the miR-30 family and have similar 
functions. Although reduced expression of members of the 
miR-30 family has been reported to be associated with breast 
tumor progression, and downregulation of miR-30a was sug- 
gested as a novel plasma marker for breast cancer and espe- 
cially for the triple negative phenotype, we did not observe a 
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difference between tumor samples and corresponding benign 
tissues in our cases [3,18,19]. 

In conclusion, our data demonstrated that miR-31, miR- 
205, miR-34a, miR-146a, miR-125b, miR-34c, and miR-181a 
were downregulated and miR-21 was upregulated more than 
1.5 -fold in TNBC tissues. In addition, our data revealed that 
lower levels of miR-205 and miR-200c have the potential to 
predict lymph node metastasis in TNBCs with high sensitivity 
and specificity. Further research with well-designed prospec- 
tive studies in larger patient groups may help to clarify the 
clinical utility of these two miRNAs as a predictive tool for 
lymph node metastasis in TNBCs. 
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